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O pewnej własności estymatorów klasycznej metody najmniejszych 
kwadratów w warunkach heteroscedastyczności składnika losowego 


O HeKOTODOM CBOJACTBE OLIEHOK KJIACCMUECKOTO METOJA HaMMEHBLIKX KBAJĄDATOB 
B yCJIOBMAX TeTepOCKEĄACTHMYHOCTM cHyUuaMHOro KOMNOHEHTA 


Von einer Eigenschaft der Astimatoren der klassischen Methode der kleinsten 
Quadrate bei der Heteroszedastizitat des Schicksalselementes 


W artykule tym zajęto się stochastycznymi właściwościami klasycz- 
nych ! estymatorów wariancji składnika losowego oraz wariancji oszaco- 
wań parametrów strukturalnych prostego modelu ekonometrycznego 
z jedną zmienną objaśniającą — w wypadku heteroscedastyczności skład- 
nika losowego. Wprowadzone zostanie pojęcie quasi-nieobciążoności esty- 
matorów, której warunki występowania poszukiwane będą każdorazowo 
przy szczególnych założeniach dotyczących modelu kształtowania się war- 
iancji składnika losowego. 


Niech będzie dany prosty model: 
(1) Y= aot aX ter  ((F1,2,..0). 


Przyjmijmy dla tego modelu klasyczne założenia: 

a) zmienna x, jest nielosowa i nie skorelowana ze składnikiem losowym, 
b) składnik losowy e+ ma wartość oczekiwaną, równą zeru oraz stałą 
wariancje o? i wszystkie autokowariancje równe 0, czyli zachodzi: 


(2a) E(e)=0 (t=1,2,...,n), 


1 Za klasyczne estymatory uważane są tu estymatory właściwe klasycznej me- 
todzie najmniejszych kwadratów. 
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o? dla t=s 

Q dla tÆs (t,8=1,2,...,n). 

Przyjmijmy na stałe, że x, jest odchyleniem od średniej z próby. Stosując 
klasyczną metodę najmniejszych kwadartów otrzymujemy estymatory 
parametrów strukturalnych modelu (1): 


(2b) E(et £s)= 


(3) ŻY 
Ao =- 
n 


X *ryt 


(4) A = 2 —— 

DAG 

(=1 o 
Przy spełnieniu założeń a) i b) estymatory (3) i (4) są na mocy twierdzenia 
Gaussa-Markowa najlepszymi nieobciążonymi estymatorami liniowymi 
parametrów m i a. 


Wariancje oszacowań (3) i (4) są równe: 
a? 
(5) D>(ap) = = 


2 
D?(1,)= RZEREM 
xi 
1 
Nieobciążonymi estymatorami wariancji składnika losowego oraz warian- 
cji (5) i (6) są odpowiednio wyrażenia: 


(6) 


M= 


(7) et 
gdzie e, są to tzw. reszty modelu, 


s? 
(8) Stay = 


2 
Sa = 


n 


2N 


f=1 


(9) 
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Należy zwrócić uwagę, że założenie b) wymaga, aby wariancja składnika 
losowego była stała, co znaczy, by zachodziło: 


D?(e1) = D?(e+) =... = D? (en). 


Przyjmijmy teraz dla modelu (1) nieco odmiene założenia: 

c) zmienna x; jest nielosowa i nie skorelowana ze składnikiem losowym, 
d) składnik e, ma wartość oczekiwaną równą 0, wariancję o, i wszystkie 
autokowariancje równe 0, czyli zachodzi: 


(10a) E(ej=0 (t=1,2,...,m), 
f 0 dla ts 
(10b) E(eies) = | 52, =ot:k, dla t=s (t,s=12,...n), 


gdzie k:>0 i Y  k=n, co oznacza, że współczynnik o? jest wariancją 
t=1 


średnią, ponieważ 


Aby założenie d) nie było tożsame z założeniem b), musi zachodzić k+7=1 
dla dwu przynajmniej różnych t. Przy tych założeniach estymatory (3) i (4) 
otrzymane za pomocą klasycznej metody najmniejszych kwadratów są 
w dalszym ciągu nieobciążone, nie posiadają natomiast własności minimal- 
nej wariancji. Dodatkową komplikacją spowodowaną przez heterosceda- 
styczność wariancji założoną w d) jest obciążenie estymatorów (7), (8), (9). 
Wariacje D?(a,) i D?fa,;), gdzie a, i a, nadal określone są przez (3) i (4), 
obecnie wynoszą: 


L m pa 
(11) GRAĆ RE: 


X. x?D? (o?) Y xłki 
(12) DJa = = nra 


D E 


natomiast ich oszacowania otrzymane za pomocą estymatorów (8) i (9) są 
obciążone. Obciążenie to dla estymatora (9) jest równe ([1] str. 310) 
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Obciążenie estymatora (8) można łatwo znaleźć korzystając z tego, że 
([1] str. 309) 


) (1—kr) e? 
E(s?)=o*| 1 + — 
(n—2))) a? 
(= 
Mamy więc 
S (1—k:) x? 
s o? | „A 
e(t 7 
(n—2) y` x? 


Odpowiednie obciążenie równa się zatem 


n 
a? >, (1—k,) x? 
= 
n(n— 2) 5 3 
RE 


Jest oczywiste, że obciążenia będą równe 0 niezależnie od wartości 
przyjmowanych przez zmienną x,, w przypadku gdyby k,=k,=k;=...= 
=kn=l, czyli gdyby zachodziły założenia a) i b). Można jednak wykazać, 
że przy założeniach c) i d) mogą występować sytuacje, w których wspom- 
niane obciążenia będą równe zeru. Przyjmując dodatkowo model, według 
którego kształtuje się wariancja składnika losowego, a więc i współczyn- 
niki k można poszukiwać takich układów wartości zmiennej objaśniającei 

` Xy dla których zachodzić będzie 


Y l=k)z8=0. 
t-—1 
Te szczególne przypadki zerowego obciążenia estymatorów (7), (8), (9) 
modelu (1) przy założeniach e) i d) oraz przyjętym modelu obiektywnej 
heteroscedastyczności nazywać będziemy dalej quasi-nieobciążonością 
albo nieobciążonością przy pewnych układach wartości zmiennej objaśnia- 
jącej. Z powyższych powodów niezbędne jest ustalenie dla dalszych roz- 
ważań modelu, według którego zmienia się wariancja składnika losowego. 
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Pawłowski [2] w przypadku estymowania funkcji popytu na dobra 
trwałego użytku przyjmuje, że wariancja składnika losowego jest funkcją 
czasu. Jak wykazują badania wielu autorów ([3] str. 256), można często 
przyjąć, że D?(e,) jest pewną funkcją bądź wartości oczekiwanej zmiennej 
objaśnianej, bądź też wartości zmiennej objaśniającej. Jeżeli na przykład 
estymujemy funkcję regresji wydatków na żywność rodzin względem 
ich dochodów, to można przyjąć proporcjonalność D? (e) do kwadratu do- 
chodu. W tej sprawie pewne ogólniejsze podejście prezentuje E. Park [3]. 
Proponuje on, aby nie zakładać apriorycznie dokładnego modelu hetero- 
scedastyczności, lecz estymować jego parametry, używając reszt otrzyma- 
nych ze wstępnego przeliczenia klasyczną metodą najmniejszych kwadra- 
tów jako oszacowania wariancji składnika losowego. Można wtedy na 
przykład estymować parametry o? i y modelu 


ae Y 
e = A; ę” 


gdzie v jest składnikiem losowym. Podejście takie — jak widać — dość 
ostrożne, unika apriorycznych założeń. Podobne sugestie wysuwa Z. Pa- 
włowski [2]. 

Wydaje się też, że możliwym do przyjęcia w pewnych sytuacjach 
założeniem jest uznanie wariancji składnika losowego za schodkową funk- 
cję indeksu t. Wiązać się to może z niejednorodnością danych zebranych 
dla celów estymacji modelu ekonometrycznego. Oczywiście bierzemy tu 
pod uwagę przypadek, gdy ta niejednorodność generuje zmienność para- 
metru stochastycznego, a nie parametrów strukturalnych estymowanego 
modelu. Należy się spodziewać, że różne modele obiektywnej heterosceda- 
styczności, które można przyjąć dla modelu (1), w zależności od jego mery- 
torycznej treści, będą implikowały różne warnuki, jakie musi spełniać 
zmienna objaśniająca, aby estymatory (7), (8), (9) można było uznać za 
quasi-nieobciążone. 

Niech będzie dalej dany model (1) przy założeniach: 

e) zmienna x jest nielosowa i nie skorelowana ze składnikiem losowym, 
f) składnik losowy et ma wartość oczekiwaną równą 0, wariancję o: i wszy- 
stkie autokowariancje równe 0, czyli: 


E(e)=0 (t= 1,2,...,0), 
Eles)= 0 dla t#s 
tês o =o?k, dla t=s (t,s=1,2,...,0), 


g) k= m't+b, gdzie mÆ+0 oraz m i b są takie, by 

23 ki=n 

1 
Wprowadzenie dodatkowego założenia g) oznacza rozpatrywanie przy- 
padku, gdy wariancja składnika losowego modelu (1) jest liniową funkcją 
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czasu. Korzystając z tego, że 
5 ki=n 
f—1 
mamy: 


" t+nb, 


1 


r 
n=m 


t 
n 
są 


1=mt +b, gdzie t =—— 


Prosta z założenia g) przechodzi więc przez punkt o współrzędnych (t, 1), 
stąd jej równanie przedstawić możemy w postaci: 

kt—1=m(t—t), 
(14) | k.=1+m(t-5. 


Esymatory (7), (8), (9) dla modelu (1) przy założeniach: e), f), g) będą 
quasi-nieobciążone, jeżeli zajdzie 


413 


(15) (1—k:)c2=0. 


~ 
m 


Wstawiając (14) do (15) otrzymujemy 


(16) m). ze(t—t)=0. 

f=1 
Równość (16) zostanie spełniona, gdy wartości zmiennej objaśniającej będą 
tworzyły szereg symetryczny. Wówczas wyrażenie 


n 

$, ze(t—t) 

1—1 
gdzie x, jest odchyleniem od odpowiedniej, będzie równe zeru. Tak więc 
w modelu (1) przy założeniach: e), f), g) warunkiem quasi-nieobciążoności 
estymatorów (7), (8), (9) jest symetryczność szeregu wartości zmiennej 
objaśniającej. 


Rozważmy teraz inny model heteroscedastyczności. 
Niech dla modelu (1) będą prawdziwe założenia: 
h) zmienna x, jest nielosowa i nie skorelowana ze składnikiem losowym, 
i) składnik losowy e+ ma wartość oczekiwaną równą 0, wariancję o? 
i wszystkie autokowariancje równe 0, czyli zachodzi: 
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E(e:)=0 (1—1,2,...,n), 
ś _ [0 dla tÆs 
(eses) = o= o?k, dla ts (t,8=1,2,...,7), 


l ki dla Ilte 
i k= |, dla tm (t=1, 2, ... m); m jest liczbą parzystą, 


oraz k;>>0, k+>0, k,g=ką, k; i k są takie aby Ð kt=n. 

t=1 
Przyjęcie założenia j) oznacza rozpatrywanie przypadku, w którym 
wariancja składnika losowego modelu (1) jest różna dla dwu równolicz- 


nych grup danych. Korzystając z tego, że > kt=n. mamy: 
f=1 


n n 
—k —k =n 
2 "Fo 2 


(17) 
ką +k,= 2. 


Estymatory (7), (8), (9), dla modelu (1) przy założeniach: h), i) j) będą 
quasi-nieobciążone, jeżeli zajdzie: 


(18) X l=kdxz2=0. 

1-1 
Korzystając z (17) oraz założenia j) przekształcamy (18) 
n ni2 n ni2 n n/2 LĄ 
Vd-k)ze=) ce], rek Y rê—k Y e?=(1—k,) (Y xe- acz) 
1-1 f—1 t=nl2+-1 11 f=n/2+-1 ET t=n/21 


Otrzymane wyrażenie będzie równe 0, jeżeli 


(=1 f-—nl2+-1 
co jest równoznaczne z równością wariancji zmiennej objaśniającej dla 
obydwu równolicznych grup danych. Aby więc estymatory (7), (8), (9) 
modelu (1) przy założeniach h), i), j) były quasi-nieobciążone, wariancje 
zmiennej objaśniającej w obydwu grupach danych powinny być sobie 
równe. 

Nakładając każdorazowo ostre założenia modelu heteroscedastyczności, 
można w pewnych przypadkach sformułować dość proste warunki quasi- 
-nieobciążoności klasycznych estymatorów wariancji składnika losowego 
oraz wariancji oszacowań parametrów strukturalnych modelu (1). Warun- 
ki te dotyczą kształtowania się zmiennej objaśniającej i są różne dla róż- 
nych typów heteroscedastyczności. Wiąże się to z tym, że wprowadzone 


9 Annales, sectio H, t. IX 
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w artykule pojęcie quasi-nieobciążoności polega na otrzymywaniu obciąże- 
nia zerowego przy pewnych układach wartości zmiennej objaśniającej. 
Jeśli ta jest zmienną ekonomiczną, układy te mogą być trudne do otrzy- 


mania. 
Należy zauważyć, że dla modelu tendencji rozwojowej postaci 


Y= aot att er (t=1,2,...,n) 


i w przypadku niejednorodności wariancji typu g) lub j) klasyczna metoda 
najmniejszych kwadratów daje zawsze quasi-nieobciążone estymatory 
wariancji składnika losowego oraz oszacowań parametrów strukturalnych. 
Zmienna czasowa spełnia bowiem obydwa znalezione przez nas warunki 
tak dla przypadku g) jak i dla j). 
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PE3IOME 


CTaTbA NOCBAINCHA CTOXACTKUECKAM NOCJIEĄCTBUAM HEOĄHOPOĄNHOCTM AKCIEPCUM 
CJIYHAKHHOTO KOMIIOHEHTA C OJHOÑ OOBACHAKIĘEJH BEJIMUKHOŃ. IIDUHATLIM METOĄROM 
ONeHKM ABJIAETCAH KIIACCMYECKMA METOĄ HaMMEHBLIAX KBAĄDATOB. Mna onperejleHuA 
HeCM€INeHHOCTM OLEHOK CTPYKTYPAJIBHBIX IAapaMeTpoB IDM HEKOTODBIX CMCTEMAX 
3HAYEHMA OOBACHAKINEŃ IIEDEMEHHOŃŃ OBIJIO BBEJĄEHO NOHATME KBAZMHECMENIEHHOCTM. 
IIOCJIeNqHAA YacTb paGOTbI NOCBAINOHA YCTAHOBJIEHMIO YCJIOBMA TOABJIEHMA JTOTO AB- 
JIeHMA NIA ONPEĄEJEHHBIX CJIYNAEB TETEpOCKENACTUHHOCTU. 


ZUSAMMENFASSUNG 


. Der Artikel befasst sich mit den stochastischen Folgen der Varianzungleichartig- 
keit des Schicksalselements des ökonometrischen Modells mit einer deutender 
Variabel. Die angenommene Methode der Ästimation ist die klassische Methode der 
kleinsten Quadrate. In der Arbeit wurde der Begriff von quasi-Unbelastetheit der 
Ästimatoren der Strukturparameter bei gewissen Wertverhältnissen von der deutender 
Variabel. Der letzte Teil befasst sich mit der Bestimmung von Bestehensbedingungen 
dieser Erscheinung für die bestimmten Fälle der Heteroszedastizität. 


